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potential P(r). In the present work, general important properties of the auxiliary 
field method are proved, such as scaling laws and independence of the results on 
the choice of P(r). The method is extended in order to find accurate analytical 
energy formulae for radial potentials of the form aP(r) + V(r), and several explicit 
examples are studied. Connections existing between the perturbation theory and 
the auxiliary field method are also discussed. 
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1. Introduction 

Auxiliary fields, also known as einbein fields, are known for a long time in quantum 
field theory. Initially they have been introduced to remove the cumbersome square 
roots appearing in relativistic theories. As an example in string field theory, let us cite 
the Nambu-Goto Lagrangian which is transformed into the Polyakov Lagrangian pQ. 
They are also of common use in many other fields of physics, such as supersymmctric 
field theories [5] and hadronic physics [3]. A particular use of the auxiliary fields is 
the transformation of a semi-relativistic kinetic energy term {\/p 2 + m 2 ) appearing 
in Salpeter type equations into an apparently non-relativistic one (p 2 /(2^)) leading 
to a simpler Schrodinger- like equation [4]. Another interesting approach involving 
auxiliary fields is the transformation of a problem containing a linear confining term 
into a new one containing a harmonic oscillator potential which can lead to analytical 
expressions [5]. 

Recently, it was realized [6] that auxiliary fields can be used as a tool to 
obtain analytical approximate expressions for the eigenvalues and the eigenstates of a 
Schrodinger equation. In this work, hereafter labelled SSB, we proposed a systematic 
method, called the auxiliary field method (AFM), which gives approximate expressions 
for the eigenenergies of a non-relativistic two-body system interacting through any 
local and central potential V(r). For special forms of this potential, an analytical 
expression is available. In SSB we studied in particular power-law and logarithmic 
potentials and proposed new energy formulae which are much more accurate than 
those found in literature up to now. 

The search for analytical solutions of the Schrodinger equation is very interesting 
and the subject of important investigations. In fact, only very few potentials give 
rise to an analytical expression for the eigenenergies valid for any values of the radial 
quantum number n and orbital quantum number I. The most famous are the harmonic 
oscillator V(r) = |wr 2 and the Coulomb potential V(r) — —n/r. 

Other potentials have analytical solutions but only for S'-waves (or in a one- 
dimensional case). This is in particular the case for the Morse potential, the Hulthen 
potential, the Hylleraas potential or the Eckart potential. There exist also potentials 
which are solvable but for particular values of their parameters. Among them, the 
Kratzer potential V(r) — 2D{\a 2 jr 2 — b/r) is famous and analytically solvable in the 
case b = a [7] . 

But even if an exact analytical expression is not available, it may be very 
interesting to have an approximate analytical expression at our disposal. In addition to 
a possible benchmark for numerical calculations, it exhibits the explicit dependence of 
the energies as a function of the various parameters and of the quantum numbers. 
Moreover, an analytical expression is always much less time consuming than the 
corresponding numerical resolution, and is thus of very great help in case of a search 
for a set of parameters for a potential relying on a chi-square best fit. 

Several methods have been invoked to find approximate analytical solutions: 
WKB method, semi-classical treatment, variational methods, perturbation theory, 
etc. In SSB we showed that the AFM is especially well suited to pursue this goal 
and presents very pleasant features. The aim of this work, which extends the results 
obtained in SSB, is essentially twofold: 

• to demonstrate a number of very general and interesting properties of the AFM, 
in particular the connections with the perturbation theory; 
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• to obtain approximate analytical expressions for a wide class of potentials, some 
of them being very important and of common use in several domains of physics. 

The paper is organized as follows. In the second section, we demonstrate a number 
of very general properties concerning the AFM. The third section is devoted to the 
application of the general theory to very specific but important potentials. The fourth 
section deals with a detailed comparison of our analytical results and the corresponding 
numerical values. Some conclusions arc drawn in the last section. 



2. General properties 

2.1. Principle of the method 

Our goal is the search for approximate analytical expressions of the eigenvalues for a 
two-body non-relativistic Hamiltonian 

H = ^ + V(r), (1) 
2m 

where m is the reduced mass, r the relative distance between particles, p the conjugate 
momentum associated with r, and V(r) any central potential. An analytical expression 
for all the eigenvalues is known explicitly only for very specific potentials P{r). The 
basic idea is to rely on such potentials. In other words, we assume that we are able 
to obtain an analytical expression e(a) for the Schrodinger equation 



it (a) 



r„2 

2m 



e(o) \a) , (2) 



in which, at this stage, a is a real parameter. 

Our method needs the introduction of an auxiliary field v, which is a priori an 
operator. We recall here the principle of the method. It consists in four steps: 

(i) We calculate the function (the prime denotes the derivative with respect to 7-) 

V'(r) 

m = (3) 

and denote by v the value of the auxiliary field which coincides with this function, 
v = K(r). (4) 

(ii) We denote by J = K _1 the inverse function of K. Thus, one has 

r = J(v). (5) 

Since both V(r) and P(r) do exhibit an analytical form, the same property holds 
for K(r). But it is by no means sure that J{v) can be expressed analytically. 
The existence of an analytical expression for J is a necessary condition for our 
method to obtain analytical expressions as a final result. 

(iii) We build a new Hamiltonian H which depends the auxiliary field v through the 
following expression 

H{v)=R{y)+g{v), (6) 

where H is defined by © and where the function g(i>) is given explicitly by 

g(u) = V(J(u)) - vP{J{v)). (7) 

This function g[y) makes a bridge between the potential P(r) for which an 
analytical expression is known and the potential V(r) for which an analytical 
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expression is a priori not known. The very important property is that v, coming 
from ([?]), cancels the variation of the new Hamiltonian, i.e. 5H {v) / 8v\ v =v = 0. 
Moreover, one has the additional crucial identity H(v) = H, as defined by ([TJ) . 
(iv) Considering now v no longer as an operator but as a pure number and taking 
into account ©, the eigenvalues of H{v) are 

E{v) = e(y)+g{y), (8) 

where e(y) are the eigenvalues of H . Then, we determine the value vq that 
minimizes E(v): dE(y)/dv\ Vo — 0. We propose to consider E{vq) as the 
approximate eigenvalues of the Hamiltonian H . In SSB, we presented a bound 
to test the accuracy of the method. In order to obtain an analytical expression 
for these eigenvalues, we must fulfill a second necessary condition: to be able to 
determine v and, then, E(i/q) in an analytical way. 

This method is completely general and a priori valid for any potential V(r). 
Nevertheless, in order to get analytical expressions we must fulfill, as we saw, two 
conditions: i) first, to be able to invert relation (j4|) in order to have access to the 
function J{v) defined by ii) second, to be able to determine vq and to calculate 
the corresponding value E(vq) in an analytical way. 

An idea for obtaining analytical expressions of the eigenenergies for an arbitrary 
potential is the following. We start with a potential P(r) — P^(r) for which the 
energies of the corresponding Hamiltonian are exactly known. We then proceed 
as above to find approximate solutions for the eigenenergies of a Hamiltonian _ff W in 
which the potential is at present V(r) = pM(r). In general, a large class of potentials 
can be treated in that way. Moreover, by comparison with accurate numerical results, 
we can even refine the expressions in order to be very close to the exact solution. 

Considering now these approximate expressions as the exact ones, we apply 
once more the AFM with P(r) = pW(r) to obtain approximate solutions for the 
eigenenergies of a Hamiltonian in which the potential is at present V(r) — P^ (r). 
Even if analytical solutions for Hamiltonian were not attainable directly with 
P{r) = pl°l(r), it may occur that they indeed are with P(r) = pM(r). Pursuing 
recursively such a procedure, one can imagine to get analytical solutions even 
for complicated potentials. Presumably, the quality of the analytical expressions 
deteriorates with the order of the recursion. 



2.2. Expression of approximate energies 

Very rare are the potentials P(r) for which an analytical solution is known for all 
radial n and orbital I quantum numbers. Among them, the harmonic oscillator (ho) 
P(r) = r 2 and the Coulomb (C) P(r) = —1/r potentials are widely used. Taking 
benefit of this opportunity, the class of power-law potentials (pi) has been studied in 
SSB. Let us consider P^(r) = sgn(A)r A with sgn(A) = A/|A| and A ^ 0. It has been 
shown in SSB that the eigenvalues of the Hamiltonian 

#A Pl) («) = |^+«sgn(A)r A (9) 
can be written under the form 

ei Pl) ( a ) = ^(a|A|) 2 /^(^) A/(A+2) , (10) 

where N\ depends on n and I quantum numbers, as well as A. This formula gives the 
exact result for the two important situations: 
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• the harmonic oscillator potential since in this case iV( n0 ) = N\ = 2 = 2n + l + 3/2; 

• the Coulomb potential for which = N\—-i = n + I + 1. 

It has been shown in SSB that, for any physical values of A (A > —2), a good form for 
N\ is given by 

N\ = b(X)n + I + c(A). (11) 

Some functions 6(A) and c(A) were proposed in SSB in order to give an approximation 
as precise as 10~ 3 for the most interesting (the lowest) values of the quantum numbers 
n and I. 

Considering the eigenenergies (jTUJ) and (fTTj) as the 'exact' ones for the 
Hamiltonian ©, one can apply the AFM to get the eigenenergies for the 
Hamiltonian (fT]). Using the recipe given in section [2TT1 with P{r) = pW(r), we find 
the following results 

E(vo) = ^ JK) A + V(J(vq)), (12) 

the optimal value i>o being determined from the equation 

N 2 

|AKJM A+2 = -±=Y x (13) 
m 

and the function J{v) coming from the relation 

\\\vJ{ v )^=V\J{v)). (14) 

Let us emphasize that J{y) depends only on the potential V(r) and not on the 
particular eigenstate we are interested in. The value vq depends both on the potential 
(through the J function) and on the state under consideration (through the N\ 
quantity). It is important to stress that the expression resulting from (|12[) suffers 
from two approximations: 

• The AFM is based on the replacement of an operator v by an optimal value uq; 
this approximation was discussed is detail in SSB and an estimation of the error 
was given; 

• Setting for N\ the value (jlip is also an approximation whose quality was discussed 
extensively in SSB. 

In the case of a harmonic oscillator A = 2 or a Coulomb potential A = — 1, only the 
first type of approximation remains. 

The application to the two solvable potentials is immediate: 

• For the Coulomb potential (A = — 1), one deduces 

£ (C) M = + V(J(^)) = ; T( \ 2 + V(J(no)). (15) 

2J(po) 2mJ((i ) 2 

The function J(fJ-) is determined by the condition 

V'(J(fi))J(fi) 2 = (i, (16) 

whereas the value /io is calculated from the transcendental equation 

fM)J(fio) = — ■ (17) 
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E^°\u ) = v I{v f + V(I(u )) = \ mI{ ^ )2 + V(I(uo)). (18) 

The function I(y) is determined by the condition 

V'{I{v))=2vI{v), (19) 
whereas the value vq is calculated from the transcendental equation 



-o/M 4 = (20) 

The comparison of the expression for the energies in both cases (fT5")) and (TT5|) clearly 
shows that there must exist a link between them. This aspect is considered in the 
next section and indeed we will prove a very interesting property. 

2.3. Switching from P^ (r) to pM(r). 

In the preceding section, we derived the approximation E x that can be obtained for 
the eigenvalues of the Hamiltonian fl} applying the AFM with the starting potential 
p( A )(r). Let us assume that, instead of the starting potential P^ (r), we are interested 
by another starting potential P«'(r). The formulae (fTI^ - lM]) apply as well, changing 
A into r\. In particular, one determines a function /(/i) depending on the auxiliary field 
fi, an optimal value /^o depending on a quantity, and the resulting energy E v (fio). 
In this last approach, let us introduce a new field v by the change of variable 

^ = Km) = MmW- a - (2i) 

Using the definition of as a function of the potential V, this new variable can be 
defined as well as 

V = v{n) = ^I{rf- x V'(I{!,)). (22) 

Using the definition of the J function HH), it is easy to show the relationship 

I(H) = JKm)). (23) 
Defining vq = v(p,o), a simple calculation shows that 

N 2 

\XWoJM X+2 = -2- = Y v . (24) 
m 

This is exactly the expression ()13[) but in which the quantity Y\ has been replaced by 
Y v - 

To achieve the demonstration, let us introduce this value ^0(1^) in the expression 
of E\(vq). Using the link between vq and [1q, it is easy to show that 

E x (u (Y v )) = ^/lo/M" + V(I(no)) = E n (v Q (Y v )). (25) 

In the expression derived from the case P^ (r) , it is sufficient to change the value N\ 
by N v to obtain the expression derived from the case P™(r). We end up with the 
very important conclusion that can be stated as a theorem: 
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If, in the expression E(N\) of the approximate energies resulting from the 
AFM with P< A )(r) ; one makes the substitution N\ — > N„ (so that E(N\) — > 
E(N V ) with the same functional form for E), one obtains the approximate 
eigenenergies resulting from the AFM with P^(r). 

In a sense, as long as we use a power-law potential ?W (r) as starting potential, 
there is a universality of the approximate AFM expression of the eigenvalue, depending 
only on the potential V(r). The only reminiscence of the particular chosen potential 
p( A )(r) is the expression of N\, as given by (jlip . This result holds whatever the form 
chosen for the potential V(r), even if we are unable to obtain analytical expressions 
for one case or the other or both. 

This property was emphasized in SSB for the particular case of a power-law 
potential V(r) — r x switching from the harmonic oscillator (A = 2) to the Coulomb 
potential (A = —1). We proved here that it is in fact totally general. It is probably 
related to a well known property in classical mechanics: one can pass from the motion 
of a harmonic oscillator to the Kepler motion by a canonical transformation. 



2.4- Scaling laws 



Scaling laws represent an important property for non-relativistic Schrodinger 
equations. They allow to give the expression for the eigenenergies (and wave functions) 
of the most general equation in terms of the corresponding eigenenergies (and wave 
functions) of a reduced equation which is much simpler to solve. 

Let us recall briefly the scaling law for the energy. Let E(m, G, a) be the 
eigenvalues of a Schrodinger equation corresponding to a system of reduced mass m 
subject to a potential of intensity G and characteristic inverse length a. The scaling 
law gives the relationship between E(m, G, a) and E{m! , G", a'). Let us start from the 
corresponding Schrodinger equations 



2m r 



GV(ar) - E(m, G, a) 



*(r) = 0, 



- — A r + G'V(a'r) 
2m 



E(m',G',a') 



*'(r) = 0. 



(26) 
(27) 



The important point is that it is the same function V(x) which appears in both 
equations. In (j27|) . let us make the change of variables r = ax and multiply it by 
X- Now, we choose the arbitrary parameters a and x m order to fulfill the conditions 
x/ijn'oi 2 ) = 1/m and era' = a. In other words, we impose the following values 

a m' / a \ 2 

<*=-, X= — [-) ■ (28) 
a m \ a' J 

With these values, (|2T|) can be recast into the form 

1 . 777/ / a. \ 2 m ' / n. \ 2 



2m 



G - 



\a'/ m \a' I \a' J 



0. 



(29) 



Equation (|26p can be recovered, provided one makes the identification G 
G' (m' /m)(a/a') 2 and a similar relation for the energies. 

The scaling law is thus expressed in its most general form as 

2 \ 



E{m > G,a) = — [- l 
m K " 



E 



m 



G— I — 

m 



(30) 
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In fact, it is always possible to define the function V(x) so that a' = 1. In what 
follows, and without loss of generality, we will apply the scaling law for energies under 
the form 

E(m,G,a) = —E(m',G'=^,l). (31) 
to \ m'a z ) 

This equality is very powerful since it is valid for the exact eigenvalues of a non- 
relativistic Schrodinger equation based on an arbitrary central potential. It allows to 
express the energy in terms of a dimensionless quantity and some dimensioned factors, 
as we will see below. 

The purpose of this section is to prove that the scaling law for energy, as expressed 
by (|3ip still holds for the approximate expressions derived from a treatment based on 
auxiliary fields. It was observed in SSB but no proof was given. The only assumption 
is that the function P is homogeneous, that is we impose the property P(ar) = a p P(r) 
which implies also P'(ar) — a p_1 P'(r). Let us denote here by M{m, G) the eigenvalues 
of the Hamiltonian 

H(G) = £-+GP(r). (32) 

In our treatment, the analytical expression for M(m, G) is supposed to be known. We 
are searching for the approximate expression of the eigenvalues of Hamiltonian ([1]) 
with V(r) — Gv(ar). We call J the inverse function of v'/P'; consequently we have 
the property v'(J(X)) = XP'{J{X)). 

Let us apply the recipe described in section 12.11 The function K(r) = 
Gav' (ar) / P' (r) — Ga p v' (ar) / P' (ar) must be identified with the auxiliary field v and 
its inversion provides r(v) with 

ar(v) = J(v/(Ga p )). (33) 

The auxiliary Hamiltonian to be considered reads 

6{y) = 7^ + vP(r) + G[v(J(jm)) - »P(J(m (34) 

with 

- (35) 

The corresponding eigenenergies are thus 

E(m, G, a; v) = M(m, v) + G[v(J(jjl)) - fiP(J(fx))]. (36) 

The next step now is the determination of the value vq which minimizes this value of 
E. A simple calculation shows that the value of /io = fi{vo) is explicitly given by 

a?M'(m, Ga> ) = P(J(j*o)), (37) 

where M'(m, v) = dM(m, v)/dv. 

Lastly, the energy E(m, G, a) we are looking for is just the value E(m, G, a; Vq). 
The final expression is therefore 

E(m, G, a) = M(m, Ga^o) + G[v{J{^ )) - fi P(J(fi ))\. (38) 



Since (|31j) is valid whatever the potential, it is in particular valid for the energies M. 
One deduces the following relations 

M (to', mGn' /(m'a 2 )) = {m/(m'a 2 ))M{m, Ga p fi' ) (39) 
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and, after differentiation, 

W (m', mGn' /{m'a 2 )) = a p M'(m, Ga P n' ). (40) 

Now let us consider the value 

E(mf, G' = (mG)/{m'a 2 ), 1) = M(m', mGfi' /(m'a 2 )) (41) 

+ {mG)l{m'a 2 )[v{J(^)) - // P(J(Mo))]- (42) 

The value of is obtained from (see (|37|) applied with a = 1) 

M'{m', mG^/{m'a 2 )) = P(J(p' )) (43) 

which, because of the property (|40|) . can be transformed into 

aPM\m,Ga p ii' Q ) = PW Q )). (44) 

A comparison between (|3"T|) and (|4"4")) immediately implies that fi' Q = /io- Thanks to 
(|39[) . this last identity inserted in (l^Tj) proves that 

E(m,G,a) = —E(m',G'=^,l). (45) 



But this result is precisely what we expect from the scaling law (see (|31|l ). 

We thus proved that scaling law for energy holds as well for the approximate 
expressions derived by the AFM. This result is almost general in the sense that it 
is valid whatever the potential under consideration, even if we are unable to obtain 
analytical solution. The only restriction is that the function P(r) is homogeneous. 
This is of no consequence since the potential P(r) = sgn(A)r A is used in practice. 

2.5. Extension of the method 

Given a potential V(r), the method for obtaining approximate solutions using auxiliary 
fields has been presented extensively in section 12.11 This method can be extended 
without difficulty for a Hamiltonian of type 

H a = + aP(r) + V(r). (46) 

One introduces an auxiliary field as before forgetting about the aP(r) contribution. 
The first 3 steps of the algorithm remain unchanged. Thus the v field is the same, 
as is the same the function g(y). The only difference arises in the expression © of 
H and H where vP{r) has to be replaced by (a + v)P{r). As a consequence, the 
corresponding energy (J8]) has to be replaced by 

E a {v) = e(a + v)+g(v). (47) 

E a (v) is an eigenvalue of Hamiltonian 

H a {v) = H{a + v)+g{v), (48) 

where H is defined by @ and where e{a + v) is an eigenvalue of Hamiltonian H(a+v). 
An eigenstate of Hamiltonians H a and H(a + v) is denoted \a + v), and we have 
e(a + v) = (a + v\H{a + v)\a + v). If vq is the value of v which minimizes (|47[) . then 
we could expect that 

E a {v Q ) = e(a + v ) + g(u ) (49) 

is a good approximation of E a , an eigenvalue of Hamiltonian J46|) . It seems that 
(|47| is very similar to (JSJ. Nevertheless, the small difference is important, because, 
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even if the determination of v from (jSJ) is technically easy, it may happen that the 
determination from (|47[) could be much more involved. 

Using the Hellmann-Feynman theorem [5] as in SSB, it can be shown that 

{a + u a \P(r)\a + v ) = P{J(v Q )). (50) 

So, J{vq) is a kind of 'average point' for the potential P(r). If functions P(r) and 
V{r) are not too different, one could expect that 

(a + v \V{r)\a + v Q ) « V (J(v )) , (51) 

and, in particular, that (see 0$ and (JSJ)) 

(a + !/ |i>|a + i/q) ~ fo- (52) 

Within this condition, the optimal value of the constant replacing the auxiliary field 
operator is close to the mean value of this operator, as mentioned in SSB. So, 
AFM could actually be considered as a 'mean field approximation' with respect to 
a particular auxiliary field which is introduced to simplify the calculations. As it is 
also shown in SSB, using (|50[) . one obtains 

E a (vo) -E a >V (J(z/ )) - (a + vo\V(r)\a + v ). (53) 

Except E a , all quantities can be computed analytically in principle. So an estimation 
of the accuracy of the approximate eigenvalue Ea^o) can be obtained. 

It is very instructive to apply (|47| for the potential P{r) — P^(r). In this 
case, the function J{y) is unchanged and still given by (fH)) (with the obvious change 
A — * rj), while the new value of vq is given, instead of (|13j) . by 

N 2 

\r 1 \{a + Vo )J{v f+ 2 = ^L (54) 

771 

and the corresponding value of the energy by 
£7(1*) = sgn ^ + ^ + ^ JM'' + V{J{v )) 



2 v ^7/(2+77) 

sgn(7?) v" I I I _ !L + K(J(i/ )). (55) 



(2 + 7 7)a + 777^0 / N; 

2 I m\rj\(a + z/ ) 



In this case, the virial theorem states that 

A + 2 

e(a + v ) = - (a + v \(a + ^ )sgn(A)r A |a + v ), (56) 

with e(z) given by (fTU|) . It is then easy to verify explicitly the equality (fBT)]) with J(vq) 
given by (|54p . These relations also imply that 

(a + f |r A |a + j/ ) = JK) A - (57) 

But, an outcome of (fSTj) is (a + fo\r\a + Uq) » J(^o) with J(^o) given by |54|) . So this 
relation is equivalent to the approximation (a + ^o|r A |a + vq) « (a+ Vo\r\a + fo) A . 

With this in mind, a potential that is worthwhile to be studied is the sum of two 
power-law potentials namely 

V(r) = sgn(n)ar n + sgn(A)&r A . (58) 

The above considerations show that the J{v) function is given by 



flvW 
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with the consequence that the optimal value is extracted from the equation 

(a + " 0) ^ ~ mh |(A +2 )/(A-,) • ( 6 °) 

In general, one does not have an analytical expression for the root of such an equation. 
This is possible only for very specific values of the powers rj and A. 

In this paper, we will study such a problem for the most favourable cases, where 
rj is chosen to give an exact expression for the eigenvalues (in practice 77 = 2 and 
r) = — 1) and where A is chosen in order to have an analytical root. 

2.6. Relation between AFM and perturbation theory 

Let us assume that V(r) = av{r) with a small enough so that av(r) <§; aP(r). If 
a is strictly 0, (|46|) and (|48|) show that v, and hence vq, vanishes. Switching on the 
potential gives a non-vanishing but small value of i>o. Indeed, v — av'(r)/P'(r), so we 
can expect that vq oc 0(a). Let us remark that J(0) can have a finite value, as it can 
be seen on (l54|) for the particular case P(r) = pW(r). 
^From results of the previous section, we can write 

E a (v) = e(a + av{J{v)) - vP(J(v)). (61) 
Using the definition the condition dE a (v)/d 

v\v=vq — implies that 
de{a + v) =e'(a + v )=P(J(v )). (62) 

It is interesting to compare this relation with (|50|) . Equation (|6T|) turns then into 

E a {vo) = e(a + u ) - u e'(a + u ) + av(J(u )). (63) 

Expanding e(a + vq) and e'(a + vq) in powers of a and keeping only the terms 0(a), 
we obtain 

E a (is )^e(a) + av(J(is Q )). (64) 

For small values of a, the contribution of av(r) can also be computed in 
perturbation. With the notations given above, we have then 

-kpert. = ( a \ H a\a) = e(a) + a(a\v(r)\a), (65) 

where \a) is an eigenstate of the Hamiltonian H(a). Since \a + v$) and |o) differ only 
by terms 0(a), we have (see ([BU]) ) 

P(J(vq)) = (a + is a \P(r)\a + v ) - (a\P(r)\a) + 0(a). (66) 

If the condition ([51) is fulfilled, we can also write 

(a\v(r)\a) n v (J(v )) + 0(a). (67) 

Using this result in (|65|) . we see that -Epert an d E a (vo) differ only by terms 0(a 2 ). 
Thus, the AFM and the perturbation theory give exactly the same results at first 
order provided V(r) do not differ too strongly from P(r). Let us remark that the 
perturbation method needs the computation of (a\v(r)\a). The AFM shows that this 
calculation can be replaced by the computation of v (J(vq)) which could be simpler in 
some particular cases. 
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3. Application to special potentials 

In SSB, we exploited the AFM for power-law potentials V(r) = br x . It was shown that 
an analytical expression for the energies exists for any value of A (we mainly focused 
our attention on values of A comprised between —1 and +2) and is given by (|10p . In 
this section, we take benefit of the remark of section l2.5l to study potentials of the form 
aP(r)^ v ' + br x . As stated above, an analytical expression is not necessarily available. 
We will examine first in which cases we do have an analytical expression, and then 
we will investigate in more details some of them. In the following, it is assumed that 
a > 0, b > and -2 < A < 2. 



3.1. Solvable potentials 

Let us examine first the case P(r) — r 2 . Application of the general method to this 
particular case leads to the following equation for the determination of v a (see ([60]) ) 

a±u =Xv£ =x , (68) 

where X = X(m,b, 7V( n °)) is some function of the parameters whose expression does 
not matter for our purpose. This equation is obviously a transcendental equation for 
which an analytical solution does not exist automatically. The only cases for which 
we are sure that an analytical solution exists is when it can be transformed into a 
polynomial of degree less than or equal to 4. In order to investigate this condition, 
let us consider 4/(2 — A) = p/q as a rational number (p and q are relatively prime). 
Calling x — Vq^ q , ([68| is recast as 

a±x q =Xx p . (69) 

All the solvable potentials should verify the conditions 1 < p < 4, 1 < <7 < 4. An 
exhaustive research of all the solvable potentials leads to the following values of the 
power A: 

A = -2, -1, -|, |, 1. (70) 
We will study in detail the cases: 

• A = — 2 because it corresponds to a centrifugal term with a real parameter instead 
of the usual 1(1 + 1) term; 

• A = — 1 because it corresponds the a simplified potential for hadronic systems 
with a short range Coulomb potential and a quadratic confinement; 

• A = 1 because this anharmonic potential is sometimes used in molecular physics. 

Now we investigate the case P(r) = —1/r. The equation corresponding to the 
determination of vq then reads for —1 < A < 2 (see (|5Tj)) ) 

(a±u )^=Z, (71) 

where again Z = Z(m,b, N^) is some unimportant function. Introducing again the 
integers p and q such that 1/(1 + A) = p/q and x — v X J q ', (171[) becomes 

(a ± x q )x p = Z. (72) 

All the possibilities to choose p and q such that p + q < 4 are suitable. The list of 
solvable potentials is given below: 
2 1 

A = -3, -2' 1 ' 2 - (73) 
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Among them we will study: 

• A = 2 because it corresponds to a potential — a/r + br 2 which is already studied 
with P(r) — r 2 . This is the only potential that can be described with either 
P(r) = r 2 or P(r) = — 1/r and this property allows very fruitful comparisons. 

• A = 1 because it corresponds to the funnel potential (Coulomb + linear) which is 
widely used in hadron spectroscopy [4]. 

It could be also interesting to introduce potentials with A < — 1 but with the 
restriction that it is repulsive at the origin (for instance, Van der Walls forces or 
Lenhard- Jones type of potentials). In this case, the equation determining vq is given 

by 



a-vo = Zv^-\ (74) 
which can be transformed with 1/(|A| — 1) = p/q and x — into 

a-x q = Zx p . (75) 
The set of all solvable potentials is provided with the list below: 

A = -5, -4, -3, --, --, -2, --, --, --, --, --. (76) 
' 2' 3' ' 4' 3' 2' 3' 4 1 ' 

In this list, we will just consider the case: 

• A = — 2 because the corresponding potential, known as the Kratzer potential, 
exhibits its spectrum under an analytical form for all values of radial quantum 
number n and orbital quantum number I. 

3. 2. Kratzer potential 

The Kratzer potential [7] is defined as 

It presents some interest as a benchmark since it is one of the rare potentials for 
which ones knows an exact analytical expression of the energies valid for any n and I 
quantum numbers. Explicitly, one has 

E(n,l) = — W - a- (78) 

n + 1/2 + + 1/2) 2 + 2ma 2 

Applying the AFM with P(r) = —1/r to this potential leads to the following 
equation giving the energies as a function of the auxiliary field v (in this section 
7V = jv(C) =n + i + \) 

2 

E ^ = -^ 2a -^~h- (79) 

The minimization with respect to v provides the value 

4ma 3 , . 

"* = wW (80) 

which, inserted in (Tf9|) . gives the desired result 

E^\n,l) = - V -2?*^ . (81) 
[2ma z + (n + I + 1) Z J 
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Let us remark that the approximate value presents the correct asymptotic 

behaviour for large n and for large I. Just to have an idea of the quality of this 
approximation, let us calculate the difference S between the denominators of E and 
It is just a matter of simple algebra to find 



6= (2n + l)(Z + l/2) 



1 



2ma? 



(Z + 1/2)* 



(82) 



Consequently, for small intensity and/or mass, ma 2 <C 1, or for large angular 
momentum, I > 1, the approximate value tends to the exact one and we have more 
explicitly 

™„2 



-(2n + l)- 



(83) 



(I + V2)' 

This behaviour is easily understandable because, under those conditions, the 
contribution due to 1/r is predominant as compared to the contribution of 1/r 2 , 
and both expressions tends towards the same exact Coulomb result. 



3. 3. Quadratic + centrifugal potential 

We consider now the potential (for an attractive centrifugal potential, not all values 
of b are relevant [9]) 

V(r)=ar 2 ±-^. (84) 

Incorporating the term ±b/r 2 into the 1(1 + l)/r 2 term already present in p 2 allows to 
get the exact eigenvalue using the same kind of arguments to those developed in the 
harmonic oscillator case [7j. Explicitly, we obtain 

E(n, l) = JJL h(2n + 1) + J (21 + l) 2 ± 8mb\ . (85) 
V 2m L 

Using the AFM with P(r) = r 2 , the energies are given by (in this section 
N = N (ho) =2n + 1 + 3/2) 



E(v) = x p^N(a + v) 1 ' 2 ± 2Vb v 1 ' 2 
V m 



(86) 



Setting Y = 2mb/N 2 , the value vq that minimizes this energy comes from a first 
degree equation and reads 



Substituting this value into the energy (|86|) . one obtains a very simple expression for 
the approximate energy 



£ (qc) (m,a,6; V) = 2y 



a(N 2 ± 2mb) 



2m 

This quantity and the corresponding exact one depend on three parameters m, 
a, b but we know that the general scaling law properties allow us to write them in a 
more pleasant form 

E(m, a, b; n, I) = \ —— e(2mb; n, I), (89) 
V 2m 



Extensions of the auxiliary field method 



15 



where e(/3; n, I) is an eigenvalue of the reduced Schrodinger equation for which the 
Hamiltonian depends now on a single dimensionless parameter f3 



H = p 2 + r 2 ±4- (90) 



r 



The exact eigenvalues of this Hamiltonian are given by 



e(/3; n, I) = 2{2n + 1) + y '(21 + l) 2 ± 4/3. (91) 

The approximate values immediately come from (|88p 

e^ c > (/?; n, I) = 2^N 2 ± /?. (92) 

One can check that the relative error between e and e(1 c ' decreases as for a 
fixed value of n and decreases as n^ 1 for a fixed value of /. This behaviour is easily 
understandable because, for large values of the quantum numbers, the contribution 
due to r 2 is predominant as compared to the contribution of 1 jr 2 and both expressions 
tends towards the same exact harmonic oscillator result. 

Let us assume that (3 <C 1; a Taylor expansion truncated to first order leads to 

e& c \(3;n,l)*2N±^. (93) 

In particular for (3 = 0, one recovers the exact value 2N, as expected. It is easy to 
check that this expression can also be obtained by perturbation theory, as expected 
from section [ 



3.4- Some reduced polynomial equations 

The potentials that remain to be studied will need the solutions of cubic and quartic 
equations. In order to simplify as much as possible the formulae, we found interesting 
to put them in a form that makes the roots as simple as possible. The corresponding 
notations will be used in the following. 

3.4-1- Cubic equation It is interesting to work with a cubic equation of the form 

x 3 + 3x - 2Y = 0. (94) 
There exists only one positive root given analytically by 



1/3 



-1/3 



F(Y)= Y+^/Y+Y 2 - Y+y/l + Y 2 . (95) 

When y « 1, one has x close to so that x 3 <C x and the behavior of the root is 
simply 

2Y 

F{Y) » — if Y < 1. (96) 

When Y ^> 1, x is large so that 3x is negligible with respect to x 3 and we have the 
following behaviour 

F(Y) w (2Y) 1/3 if Y > 1. (97) 
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3.4-2. Quartic equation The quartic equation which gives the most pleasant form for 
the roots is 

4a; 4 ± 8a; - 3Y = 0. (98) 
For each sign, there exists only one positive root given analytically by 

G±(Y) = T l -^VW) + \^{V{Y))-^ V (Y), (99) 

with 

v(Y) = (2 + y/l + y 3 ) 1/3 - y (2 + ^4 + r 3 ) 1/3 . (100) 

When Y -C 1, one has x A <C x and the behaviour of the roots is simply 

G + {Y)^— if y«i, (101) 

and 

G_(y)»2 1 / 3 + — if y<i. (102) 

8 

When Y ^> 1, a; is large so that 8a; is negligible with respect to 4a; 4 and we have the 
following behaviour 

G ± (Y) w (3Y/4) 1 / 4 if y>l. (103) 

5.5. Anharmonic potential 

The potential under consideration reads 

V(r) = ar 2 + 2br. (104) 

The interest for such a potential is discussed in section l3~Tl Obviously, one must take 
P(r) = r 2 in the AFM and the energies to be considered are given by (in this section 
N = 7V( h °) =2n + / + 3/2) 

nr b 2 

EM = J— N{a + v) 1/2 + — . (105) 

V TO V 



Let us introduce the parameter 

8 ( N 2 
3° lm6 4 



and the new variable 

*-fj K < 107 > 

The equation that leads to the minimization of the energy is then 

4a; 4 - 8a; - ZY = 0. (108) 

This is the reduced quartic equation studied in the preceding section. The solution is 
given by Xq(Y) = G_(Y) (see V Substituting this value into the energy (| 1 05|) and 
making a little algebra leads to the desired approximate energy 

E^(m,a,b;n,l) = ^Y (g 2 _{Y) + ) , (109) 

with Y given by (fT06|) . 
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As in the previous case, this quantity and the corresponding exact one depend on 
three parameters m, a, b but the general scaling law allows us to write them in terms 
of a reduced quantity depending on a single parameter [3 

[2a ( 3b 2 [3m \ 
E(m,a,b;n,l) = — e ^; n, I j , (110) 

where e(/3; n, I) is an eigenvalue of the reduced Schrodinger equation for the 
Hamiltonian 

r> 2 

H=^j + 3r 2 + 8yff3r. (Ill) 
The approximate value corresponding to this reduced equation follows from (|109[) 

e ( an )(/3;n,/) = 2/3r(^G 2 _(y) + ^y), Y = (jj ^ . (112) 

The parameter (3 could also be associated with the quadratic potential, but this case 
less interesting is not considered here. 

Let us assume that /3 <C 1; a Taylor expansion truncated to first order leads to 



£ (an) 



(P;n,l)^V3N + 4 ] J^. (113) 

In particular for (3 = 0, one recovers the exact value s/3N, as it should be. This result 
comes also from the perturbation theory. 

The limit /3 — » oo is not physically relevant, but it is interesting to consider it in 
order to check the formula. In this limit, we find 

e (an) (/3;n,0 = 3(4/37V 2 ) 1/3 + O (/T 1/3 ) . (114) 

The dominant term is the result expected for a pure linear potential, as given by (|10p . 

3.6. Quadratic + Coulomb potential 

In this section, we study the quadratic + Coulomb potential defined as 

V(r) = ar 2 - -. (115) 

r 



The interest for such a potential is discussed in section I3TT1 Let us illustrate the AFM 
method with the option P(r) = r 2 so that N = N^ ) =2n + l + 3/2. The energies 
depending on the auxiliary field v are given in this case by 

2~ , ,i/o „/& 2 ^ xV3 



E(v) = x -N(a + v)V*-3 — . (116) 
m V 4 



Let us introduce the parameter 

8A^ 2 [Aa^ 1/3 

and the new variable 

( , 

v ) 



2a \ 1/3 

) • (H8) 
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The equation that leads to the minimization of the energy is then 



4x 4 



ix - 3Y = 0. 



(119) 

This reduced quartic equation was studied previously. The solution is given by 
xq(Y) = G+(Y) (see (|99|) ). Substituting this value into the energy (|116p . one is 
led, after some manipulations, to the desired approximate energy 

U2\ 1/3 



E 



(qC), 



3 fab' 



m, a, b;n,l) = — 



Y 



Gi(Y) G+(Y) 



(120) 



with Y given by (|117j) . 

One can check that, starting with P(r) = — 1/r and performing the same kind of 
algebra, the approximate energy is given by exactly the same equation as (|120[) . but 

this time with the Y parameter given by dTTTJ) in which 7V( ho ) is replaced by 

This property is a very nice check of the general prescription demonstrated in section 

EH 

The quantity (|120p and the corresponding exact one depend again on three 
parameters m, a, b but the general scaling law allows us to write them in terms of a 
reduced quantity depending on a single dimensionless parameter (3. One can imagine 
two formulations depending on whether (3 is part of the quadratic contribution or of 
the Coulomb contribution: 



E(m, a, b; n, I) = 4-t / - — e — 



2a / 1 {Mm 3 b 4 



3 m 

E(m,a,b;n,l) = —— r] 
lb 



4 V 



1/6 



;n,l 



x 1/6 



54m 3 6 4 



n, / 



(121) 



(122) 



The e and 77 energies are the eigenvalues of the reduced Schrodinger equations for the 
respective Hamiltonians H t et H v : 



3p^ 
16 

3p 2 V2 



^3/2 



(123) 
(124) 



16 r 

The approximate values corresponding to these reduced Hamiltonians follow from 



e (qC) (/W) = M 



^ C )(/3>,0 = 



3/3'^ 



G^(F) G+(F) 
Y 4 



G 2 + (y) G+(y) 



Y = (N/3'f 



(125) 
(126) 



Let us assume that, (3 <C 1, that is V(r) <C P(r): 

A Taylor expansion truncated to first order for the formulation based on the e 
form leads to 



,(qC) 



I 2/3 3 



(127) 



In particular for (3 = 0, one recovers the exact value ViN/A if JV = 7V( h °). 
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• The Taylor expansion for the formulation based on the r\ form gives 

In particular for 0' = 0, one recovers the exact value — 8/(3iV 2 ) if N — . 

In all cases, the approximate formulae resulting from AFM agree with the result of 
perturbation theory. 

The limit (3 — > oo is interesting to consider in order to check the formulae: 

e ^\f3;n,l) = -^ + 0(f3- 2 )- (129) 
If N = N<- C \ the dominant term is the exact result for a Coulomb potential. 

^ C )(/3';n,0 = ^/3' 3 7V + o(/3' 3/2 ). (130) 

If TV = 7V( n0 ) , the dominant term is the exact result for a quadratic potential. 

Let us emphasize the point that both (| 123[) and (|124[) can be related with the 
scaling laws developed in section [2T4l As a consequence, it can be shown that both 
the exact eigenvalues and the AFM approximate ones fulfill the relation 

e(f3;n,l) = ^r](l/ f3;n,l). (131) 
3. 7. Funnel potential 

In this section, we are concerned with the funnel potential defined as 

V(r) = ar--. (132) 
r 

This potential is particularly important and its interest is discussed in section 13.11 
Let us recall that it is widely used in hadronic spectroscopy and corresponds to 
a linear confinement coupled to a short range Coulomb contribution 4J. Finding 
approximate analytical values for the energies corresponding to this potential is thus 
a very interesting question. To our knowledge such formulae are not proposed in 
literature. 

Naturally, one must apply the AFM method with P(r) = —1/r so that N — 

N (C) 

= n + / + 1. The energies depending on the auxiliary field v can be calculated 
following the prescription detailed in section f2. 51 Explicitly, one finds 



2Af 2 

Let us introduce the parameter 



E{y) = - ^±^ + 2V-a^. (133) 



3 2 / 3a 



(135) 



and the new variable 

The equation that must be solved to minimize the energy is 

x 3 + 3x - 2Y = 0. (136) 
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This reduced cubic equation was studied previously. The solution is given by 
Xo(Y) = F{Y) (see (IMI) ). Inserting this value into the energy (|133|) . one finds the 
expression of the approximate energy 



E^ (to, a, b; n, I) = V Sab 



Y 



(137) 



_F 2 (Y) F(Y)_ 

with Y given by (|134[) . Let us mention that another form of this equation can be 
found thanks to the following relation 

Y 2 1 

„, r — = sinhfl , (138) 

F 2 {Y) F(Y) 4sinh0 V ' 

with the change of variables Y = sinh(30). 

The quantity (|137p and the corresponding exact one depend again on three 
parameters to, a, b but the general scaling law allows us to write them in terms 
of a reduced quantity depending on a single dimensionless parameter (3. In hadronic 
physics, the dominant interaction between a quark and an antiquark is a confining 
linear potential [4 J. So, although the linear potential has no analytical exact solution 
for all values of the quantum numbers, it is also interesting to consider two formulations 
depending on whether (3 is part of the linear contribution or of the Coulomb 
contribution: 



E(m, a, b; n, I) = 3 



2m 



1/3 



'4m 2 o 3 



1/4 



\ V 27a 

r?< h n 2mb2 (( 27a V /4 / 
E(m, a, b;n,l) = _ r? I I I ;n,l 



;n,l 



(139) 



(140) 



The e and rj energies are the eigenvalues of the reduced Schrodinger equations for the 
respective Hamiltonians H e et H n : 

,2 „ ^4/3 



3 3 



H v -- 



— +P' A r. 
r 



(141) 



(142) 



The approximate values corresponding to these reduced Hamiltonian follow from (|137|) : 



Y 



F 2 (Y) F{Y) 



Y 



Y = 



Y = (N/3'Y 



-,%>,<> = 3*/ V ny) f(y) 

Let us assume that /3< 1, that is V(r) <C P(r): 
• The Taylor expansion for the formulation based on the e form gives 

N 2 / 3 



/ /3 4 \ 1/3 

2 2/3 UiV 2 J ' 



(143) 
(144) 

(145) 



In particular for (3 = 0, one recovers the value expected for a pure linear potential, 
as given by (fi"0|) . 

The Taylor expansion for the formulation based on the rj form gives 

35/3 2N 2 /3' 4 



/)(/?'; n,l) 



AN 2 3 4 / 3 ' 
In particular for (3' = 0, one recovers the exact value 



(146) 



-3 5 / 3 /(47V 2 
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In all cases, the approximate formulae resulting from AFM agree with the result of 
perturbation theory. 

The limit (3 — > oo is interesting to consider in order to check the formulae: 

e( f )(/3;n,0 = -|^ + o(r 4/3 ). (147) 
The dominant term is the exact result for a Coulomb potential. 

/)(/3';n,0= Q/3' 4 iv) 2/3 + 0(/3' 4/3 ). (148) 
The dominant term is the result expected for a pure linear potential, as given by 

ma. 



Here again, both (|141[) and (|142[) can be related with the scaling laws. As a 
consequence, it can be shown that both exact eigenvalues and AFM approximate ones 
fulfill the relation 

£8/3 

e(/W) = ^a/s'KV/W)- (149) 



4. Comparison with numerical results 

In the previous section, approximate analytical forms for eigenvalues of several 
Hamiltonians were found. The formulae depend on the quantum numbers n and I 
through a factor N. This number could be taken as N^°\ ) or even N\ if 
the potential P(r) chosen is the power-law potential P^(r) (see section l2~3j) . If we 
look, for instance, at the Hamiltonian (|123p . it is clear that it reduces to a harmonic 
oscillator when (3 = 0. In this case, the choice N = N^ 10 ^ gives the exact result. 
When [3 — > oo, the Coulomb part dominates and the choice N = 

N (C) 

is expected to 

yield the exact result asymptotically. 

All dimcnsionlcss Hamiltonians considered above depend on a parameter (3. The 
variation of the eigenvalues being smooth for the variation of /?, we can assume that 
the number N, giving the optimal values for a selected set of these eigenvalues, is also 
a smooth function of (3. From considerations above, the functional form 

N{(3) = b((3)n + I + c{(3) (150) 

seems reasonable. If, in the limits — ► and (3 — > oo, the Hamiltonian reduces to a 
known form, the values N(0) and iV(oo) can be computed. But for finite values of j3, 
we cannot predict the correct behaviour. It is then necessary to focus our attention 
on numerical solutions. 

Very accurate eigenvalues enum(/3; n, I) for Hamiltonians defined above can be 
obtained numerically with the Lagrange mesh method [lOj . It is very accurate and 
easy to implement. In order to find the best possible values for coefficients d(j3) (d 
stands for & or c), we will use the following measure 

1 3 3 

M = ^EE( £ num(/3 ;fl ,l) - e a pp(/?; n, i)f , (151) 

n=0 1=0 

where eapp(/3; n, £) are values obtained from our approximate formulae. Other choices 
are possible but we find this one very convenient. The analytical form eapp depends 
on N(f3) which depends on coefficients d. For each value of (3, optimal values for 
the d coefficients, d m ^ n (f3), can be determined by minimizing x(/3). Then, with a 
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set {e? mm (/?)} for a given set {/?}, a functional form d^{(3) can be fitted with the 
following measure 

x(d)=J2(d min ([3)-d m (P)) 2 . (152) 

Again, other choices are possible but we find this one very convenient. We will now 
try to determine the best form of coefficients d(/3) for some of the potentials studied 
above. 

4--1- Improvement for anharmonic potential 

The dimensionless Hamiltonian proposed here for the anharmonic potential is 

H=^ + 3r 2 + 8y^r. (153) 

Approximate eigenvalues are given by (Tm)) . With ([150]) , the exact result will be 
obtained for 6(0) = 2 and c(0) = 3/2. When — ► oo, the linear part dominates 
and, from results obtained in SSB, we could expect that 6(oo) ~ tt/VS ~ 1.814 and 
c(oo) w \/37r/4 w 1.360. By minimizing our measure x{P)i we found the optimal 
values of b((3) and c(j3) for several finite values of f3. The results are plotted with dots 
on figure [1] One can clearly see the smooth transition between two domains for zero 
and infinite values of (3. This corresponds to the transition between two potentials 
not too different, a linear one and a quadratic one. 



b([3) 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 

Figure 1. Best values of the coefficients 6(/3) and c(/3) to parametrize the 
eigenvalues of Hamiltonian (1 1 53 ft : numerical fit with (1151 ft [dots]; functions 11154ft 
with set 1 of parameters from table [T] [solid line] ; idem with set 2 [dashed line] . 



We tried to fit the numerical points with various functions and found that 
the best result is obtained for sections of hyperbola, both for 6(/3) and c(/3). The 
parameterization retained is 

5 /9 = ~1TZ ' c ^ = ' ( 154 ) 

P+P3 P + Q3 

With this choice, 6(0) = P2/P3, 6(00) = pi, c(0) = 92/93 and c(oo) = q\. Two different 
fits are presented in table [1] For set 1, only one parameter is free for each coefficient 
and the following constraints are imposed: 6(0) = 2, 6(00) = c(0) = 3/2, 

c(oo) = s/3w/4. For set 2, all parameters are free; the results found by minimization 
give: 6(0) « 1.990 close to 2, 6(00) = 1.826 close to 1.814 (tt/\/3), c(0) w 1.496 close 
to 1.5, c(oo) = 1.381 close to 1.360 (V3tt/4). 
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Table 1. Values of parameters pi and qi for 1 1154 l l. Fixed parameters are marked 
by a *. 





Pi 


P2 


P3 


51 


52 


<?3 


Set 1 
Set 2 


1.826 


2 P 3* 

1.485 


0.835 
0.747 


1.381 


3g 3 /2* 
0.333 


0.445 
0.222 



The quality of the fits can be appraised by examining the values of x(P) shown in 
table [2J It is clear that allowing a /3-dependence for the coefficients b and c improves 
greatly the approximate eigenvalues. Let us remark that the fit with three parameters 
is only slightly better than the fit with only one parameter, and that x(P) increases 
with [3 for the choice N = N^°\ as expected since the potential deviates more and 
more from a pure quadratic one. 



Table 2. Values of x(P) f° r eigenvalues of Hamiltonian (1 1 53 ft as a function of /3, 
for various parameterization of N = b(/3)n + Z + c(/3) (see table [TJ. Af( no ' and set 
1 give the exact result for /3 = 0. 






AT(ho) 


Set 1 


Set 2 


0.1 


5.9 10~ 2 


6.2 10~ 3 


4.0 10 -3 


1 


0.46 


2.7 10~ 3 


1.9 10~ 3 


10 


2.8 


1.8 10" 2 


5.4 10" 3 



4-2. Improvement for quadratic + Coulomb potential 

To study the case of the quadratic + Coulomb potential, we choose the following 
dimensionless Hamiltonian 

H44-^. (155) 
16 4 r 

Approximate eigenvalues are given by (|125p . The exact result will be obtained for 
6(0) = 2 and c(0) = 3/2. When (3 — > oo, the Coulomb part dominates and we could 
expect that 6(oo) = 1 and c(oo) = 1. By minimizing our measure x(fl)-> we found 
the optimal values of b((3) and c(f3) for several finite values of (3. The results are 
plotted with dots on figure [21 One can clearly see the smooth transition between the 
harmonic oscillator region near zero and the asymptotic Coulomb region. In this case, 
the transition occurs between two very different potentials: r 2 and — 1/r. It is clear 
that the asymptotic region is more rapidly reached for c(f3) than for b{0). When the 
quantum numbers n and I increase, the size of the eigenfunctions grows and the states 
become less sensitive to the Coulomb part. So, we can understand that the influence 
of [3 on the coefficient b((3) is less significant than on c((3). 

We tried to fit the numerical points with various functions and found that the 
best result is obtained for exponential functions with a cubic argument, both for b{(3) 
and c{[3). The parameterization retained is 

b(f3) = l+ Pl cxp(- P2 (f3- P3 ) 3 ), c(i3) = l + qi cx P (-q 2 (p~q 3 ) 3 ). (156) 

With this choice, 6(oo) = 1 and c(oo) = 1. Two different fits are presented in table [3] 
For set 1, only one parameter is free for each coefficient and the following constraints 
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0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 

Figure 2. Best values of the coefficients 6(/3) and c(/3) to parametrize the 
eigenvalues of Hamiltonian l|155| l: numerical fit with (|151|l [dots]; functions l|156jl 
with set 1 of parameters from table \3\ [solid line] ; idem with set 2 [dashed line] . 



are imposed: 6(0) = 2 and c(0) = 3/2. For set 2, all parameters are free; the results 
found by minimization give: 6(0) « 1.991 close to 2 and c(0) « 1.495 close to 1.5. 

Table 3. Values of parameters pi and qi for l|156j l. Fixed parameters are marked 
by a *. 





Pi 


P2 


P3 


91 


92 


93 


Set 1 
Set 2 


1* 
0.990 


0.093 
0.119 


0* 
0.161 


1/2* 
0.496 


2.414 
1.373 


0* 
-0.136 



The quality of the fits can be appraised by examining the values of x(P) shown in 
table |U Again, allowing a /3-dependence for the coefficients 6 and c improves greatly 
the approximate eigenvalues. The fits, with one parameter and three parameters, 
give nearly the same result. x(P) increases with f3 for the choice N = N^°\ as 
expected since the potential deviates more and more from a pure quadratic one. On 
the contrary, x(P) decreases with (3 for the choice N = as expected since the 

potential comes closer to a pure Coulomb one. The values of x{P) are large because 
the considered values of (3 are not in the asymptotic region of b(/3), as one can see on 
figure [21 It has been checked that x(P) ~ ¥ also for sets 1 and 2 for large values of (3. 

Table 4. Values of x{P) f° r eigenvalues of Hamiltonian d 1 55 ft as a function of fj, 
for various parameterizations of N = b({3)n + / + c(/3) (see table J\K ) gives 
the exact result for /3 = 0, 

N (C) 

gives the exact result for j3 = do, set 1 gives 
the exact result for both j3 = and /3 = oo, and set 2 gives the exact result for 



p = 


oo. 











Ar(ho) 


N (C) 


Set 1 


Set 2 


0.5 


9.6 10~ 3 


0.90 


2.6 10~ 3 


2.3 10" 3 


1 


0.11 


0.75 


1.9 10" 2 


1.8 10" 2 


2 


3.1 


0.45 


0.12 


0.12 
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4-3. Improvement for funnel potential 



For the case of the funnel potential, we choose the following dimcnsionless Hamiltonian 



3 3 



/3 4 / 3 



(157) 



Approximate eigenvalues are given by (|143p . This choice is motived by three reasons: 

• For hadronic problems [3], physical values of (3 varies from to about 1.5; 

• A comparison is possible between Hamiltonians (|157p and (|155[) which differ by 
the 'confinement' part. 

• Contrary to previous cases, (3 = does not correspond to a full analytical case; 
For [3 = 0, the potential is pure linear one. So, from the study performed in SSB, we 
can expect that 6(0) = 7r/v3 » 1.814 and c(0) = y/Sn/A w 1.360 is a good choice. 
When (3 — > oo, the Coulomb part dominates and we could expect that 6(oo) = 1 and 
c(oo) = 1. By minimizing our measure x(P)i we found the optimal values of b(f3) and 
c(f3) for several finite values of (3. The results are plotted with dots on figure [3] One 
can see a smooth transition similar to the previous case. Again, for probably the same 
reason, the asymptotic region is more rapidly reached for c{(3) than for b(f3). 



1.8 




c(/3) 



1.3 



1.2 



1.1 



0.5 1 1.5 2 2.5 3 3.5 



Figure 3. Best values of the coefficients 6(/3) and c(/3) to parametrize the 
eigenvalues of Hamiltonian JT57J: numerical fit with (1151 ft [dots]; functions II158H 
with set 1 of parameters from tabic \5\ [solid line] ; idem with set 2 [dashed line] . 



We tried to fit the numerical points with various functions and found that the best 
result is obtained for a Gaussian function both for b{(3) and c([3). The parameterization 
retained is 

b(f3) = l+p 1 cxp(-p 2 2 (f3-p 3 ) 2 ), c({3) = l + q 1 cxp(-q 2 2 (P-q 3 ) 2 ). (158) 

Let us note the difference with the previous case, for which the argument of the 
exponential was cubic in (3 and not quadratic. With this choice, 6(00) = 1 and 
c(oo) = 1. Two different fits are presented in table O For set 1, only one parameter is 
free for each coefficient and the following constraints are imposed: 6(0) = and 
c(0) = V3V/4. For set 2, all parameters are free; the results found by minimization 
give: 6(0) » 1.774 close to tt/VS (1.814) and c(0) w 1.367 close to VSn/A (1.360). 

The quality of the fits can be appraised by examining the values of xift) shown 
in table [6] and the values of approximate results compared with exact ones presented 
in table [7J As in the two previous cases, allowing a /^-dependence for the coefficients 
6 and c improves greatly the approximate eigenvalues. The fits with one parameter 
and three parameters give also nearly the same results and the behaviour of x(P) f° r 
this case and the previous one are very similar. 
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Table 5. Values of parameters pi and qi for H158I I, Fixed parameters are marked 
by a *. 





Pi 


P2 


P3 


9i 


32 


93 


Set 1 
Set 2 


7T 1 * 

0.783 


0.416 
0.459 


0* 
0.237 


"\/37r -1 * 

4 1 

0.369 


1.245 
1.168 


0* 
-0.062 



Table 6. Values of x(P) f° r eigenvalues of Hamiltonian JT57J as a function of 0, 
for various parameterization of N = b(/3)n + I + c(/3) (see table IB1 . 

AT(C) 

set 1 

and set 3 give the exact result for /3 = oo. 



p 


N (C) 


Set 1 


Set 2 


0.5 


0.17 


4.2 10" 4 


3.6 10" 4 


1 


0.15 


2.9 10" 3 


2.9 10" 3 


2 


0.11 


1.7 10~ 2 


1.7 10~ 2 



Table 7. Eigenvalues e(/3o,n,l) of Hamiltonian J 1 571 ) with /?o = 0.5, for some 
sets (n, I). First line: enum(0o',n, £) from numerical integration; second line: 
e^((3 ;n,l) given by Jl43l with N(/3) denned by set 1; third line: e& (/3 ;n,l) 
given by l(T43) l with N(J3) = N^ C \ 



1 


e(Po;a,l) 


e(A>;i.Q 


6(A>;2,z) 







0.39711 


1.11714 


1.64558 


2.09628 




0.42779 


1.16223 


1.68099 


2.12205 




0.26827 


0.79105 


1.15440 


1.45987 


1 


0.90598 


1.45955 


1.92580 


2.34167 




0.88794 


1.46673 


1.93564 


2.34911 




0.79105 


1.15440 


1.45987 


1.73269 


2 


1.25749 


1.74247 


2.17133 


2.56288 




1.23307 


1.73892 


2.17323 


2.56506 




1.1544 


1.45987 


1.73269 


1.98358 


3 


1.55457 


1.99727 


2.39917 


2.77168 




1.52908 


1.98937 


2.39764 


2.77183 




1.45987 


1.73269 


1.98358 


2.21833 



4- 4- General considerations 

^From the results of this section, it is clear that a good choice for the function N((3) 
can greatly improve the accuracy of the energy formulae. Unfortunately, the best 
functional form cannot be theoretically predicted. In the cases studied in this paper, 
the same form can be given to both coefficients b(f3) and c(/3) for a given Hamiltonian, 
but with different parameters. The behaviour for j3 — and (3 — > oo can sometimes 
be exactly computed. 

Obviously, the parameters for coefficients b((3) and c([3) depend on the points 
used for the fit but also on the particular choice of functions (| 15 1|) and (|152|) . Other 
definitions-relative error instead of absolute error, different chosen quantum numbers 
or different chosen values of (3 in summations-would have given other numbers slightly 
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different. It is worth noting that, for each case studied, values of xW) obtained directly 
with coefficients d mm (/3) and d^AfS) for the set 2 are generally very close. 

5. Conclusions 

The AFM was proposed in SSB as a tool to compute approximate analytical solutions 
of the Schrodinger equation and then applied to the case of power-law radial potentials. 
The basic idea underlying this method is to replace an arbitrary potential V(r), for 
which no analytical spectrum is known, by an expression of the type v P(r) + g{v), 
P(r) being a potential for which analytical cigcncncrgies can be found, v the auxiliary 
field and g(y) a well-defined function of this extra parameter. This auxiliary field is 
such that its elimination as an operator leads to the original Hamiltonian. If v is seen 
as a number however, analytical eigenenergies and eigenstates can be found, and the 
auxiliary field is eventually eliminated by a minimization on the eigenenergies. The 
approximation in the AFM comes thus from the use of the auxiliary field as a number 
rather than an operator. 

In the present work, we have further investigated the AFM and obtained results 
that fall mainly in two categories: general properties of the AFM and analytical 
resolution of the Schrodinger equation with potentials of the form ar x ± br v . Let us 
first summarize the general properties of the AFM that we have proved: 

• The analytical expressions that are obtained for the eigenenergies by using 
the AFM preserve the general scaling laws ([30|) in the case where P(r) is 
homogeneous. Thanks to this feature, the number of relevant parameters in the 
Hamiltonian can be seriously reduced in order to simplify the problem. 

• Let Pi (r) and P2 (r) be two power-law potentials whose eigenenergies respectively 
read Ei(Ni) and £'2(^2), where N± and N% are terms containing the radial and 
orbital quantum numbers. Then, if we apply the AFM to find the eigenenergies 
of potential V(r) with both Pi(r) and P2 (r), the final results will have the same 
functional form, depending either on N\ or on N2 following the case. 

• Perturbation theory can be reformulated within the AFM. If a is a small 
parameter and if V(r) — o~v(r) is a potential that can be treated in perturbation, 
then (v(r)) can be equivalently replaced by v(J(vq)) with J(Vo) the average point 
defined by (|50|) . at the first order in a. Such a property could be useful in some 
cases since no integration is needed with the AFM. 

Finally, it is worth summing up the results that we obtained by solving the 
Schrodinger equation with potentials of the form ar x ± br 11 : 

• Analytical formulae have been found for several potentials that are relevant in 
various domain of physics: anharmonic, quadratic plus Coulomb and funnel. 
To our knowledge, it is the first time that such formulae are found. A 
comparison with the Kratzer and the quadratic plus centrifugal potentials which 
are analytically solvable is also performed. 

• By using the scaling laws, it appears that only one dimensionless parameter, 
denoted as f3 in this work, 'controls' the features of the various spectra. The 
approximate formulae for eigenenergies can be a complicated function of (5 and 
of a number N containing the radial and orbital quantum numbers. This number 
is determined by the potential P(r) chosen and is a priori independent of j3. 
Nevertheless, a drastic improvement of the formulae we computed can be obtained 
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by replacing this term TV by a function of the form b((3)n+£+c(f3), where b{0) and 
c{(3) have to be fitted on the eigenenergies coming from a numerical resolution of 
the Schrodinger equation. This has finally led us to very accurate formulae for 
all the various potentials we studied. 

As an outlook, we mention that more complicated potentials could be studied 
with the AFM. In particular, we plan to apply this method to potentials involving an 
exponential, like the Yukawa potential for example. Such a work is in progress. 
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